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Abstract. We consider the problem of assigning radii to a given set of 
points in the plane, such that the resulting set of circles is connected, and 
the sum of radii is minimized. We show that the problem is polynomially 
solvable if a connectivity tree is given. If the connectivity tree is unknown, 
the problem is NP-hard if there are upper bounds on the radii and open 
otherwise. We give approximation guarantees for a variety of polynomial- 
time algorithms, describe upper and lower bounds (which are matching 
in some of the cases), provide polynomial-time approximation schemes, 
and conclude with experimental results and open problems. 

Keywords: intersection graphs, connectivity problems, NP-hardness prob- 
lems, approximation, upper and lower bounds 

1 Introduction 

We consider a natural geometric connectivity problem, arising from assigning 
ranges to a set of center points. In a general graph setting, we are given a weighted 
graph G = (V, E). Each vertex v e V in the graph is assigned a radius r v , and 
two vertices v and w are connected by an edge in the connectivity graph 
H = (V, F), if the shortest-path distance d(v, w) in G does not exceed the sum 
r v +r w of their assigned radii. In a geometric setting, V is given as a set of points 
P = {pi, . . . ,p n } in the plane, and the respective radii rj correspond to circular 
ranges: two points pi, pj have an edge in the connectivity graph, if their circles 
intersect. The Connected Range Assignment Problem (CRA) requires an 
assignment of radii to P, such that the objective function R — J2i , a = 1 is 
minimized, subject to the constraint that H is connected. 

Problems of this type have been considered before and have natural motiva- 
tions from fields including networks, robotics, and data analysis, where ranges 
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have to be assigned to a set of devices, and the total cost is given by an objective 
function that considers the sum of the radii of circles to some exponent a. The 
cases a = 2 or 3 correspond to minimizing the overall power; an example for 
the case a = 1 arises from scanning the corresponding ranges with a minimum 
required angular resolution, so that the scan time for each circle corresponds to 
its perimeter, and thus radius. 

In the context of clustering, Doddi et al. [f , Charikar and Panigraphy [5], 
and Gibson et al. [9] consider the following problems. Given a set P of n points 
in a metric space, metric d(i, j) and an integer k, partition P into a set of at most 
k clusters with minimum sum of a) cluster diameters, b) cluster radii. Thus, the 
most significant difference to our problem is the lack of a connectivity constraint. 
Doddi et al. [7] provide approximation results for a). They present a polynomial- 
time algorithm, which returns O(k) clusters that are 0(log(^))-approximate. For 
a fixed k, transforming an instance into a min-cost set-cover problem instance 
yields a polynomial-time 2-approximation. They also show that the existence 
of a (2 — e)-approximation would imply P — NP. In addition, they prove that 
the problem in weighted graphs without triangle inequality cannot be efficiently 
approximated within any factor, unless P — NP. Note that every solution to 
b) is a 2-approximation for a). Thus, the approximation results can be applied 
to case a) as well. A greedy logarithmic approximation and a primal-dual based 
constant factor approximation for minimum sum of cluster radii is provided by 
Charikar and Panigraphy [5] . In a more geometric setting, Bilo et al. [5] provide 
approximation schemes for clustering problems. 

Alt et al. PQ consider the closely related problem of selecting circle centers 
and radii such that a given set of points in the plane are covered by the circles. 
Like our work, they focus on minimizing an objective function based on rf 
and produce results specific to various values of a. The minimum sum of radii 
circle coverage problem (with a = 1) is also considered by Lev-Tov and Peleg 
[10] in the context of radio networks. Again, connectivity is not a requirement. 

The work of Clementi et al. [B] focuses on connectivity. It considers minimal 
assignments of transmission power to devices in a wireless network such that 
the network stays connected. In that context, the objective function typically 
considers an a > 1 based on models of radio wave propagation. Furthermore, 
in the type of problem considered by Clementi et al. the connectivity graph 
is directed; i.e. the power assigned to a specific device affects its transmission 
range, but not its reception range. This is in contrast to our work in which we 
consider an undirected connectivity graph. See [Sj for a collection of hardness 
results of different (directed) communication graphs. 

Carmi et al. [4] prove that an Euclidean minimum spanning tree is a constant- 
factor approximation for a variety of problems including the Minimum-Area 
Connected Disk Graph problem, which equals our problem with the different ob- 
jective of minimizing the area of the union of disks, while we consider minimizing 
the sum of the radii (or perimeters) of all circles. 

In this paper we present a variety of algorithmic aspects of the problem. In 
Section 2 we show that for a given connectivity tree, an optimal solution can be 
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computed efficiently. Section 3 sketches a proof of NP-hardness for the problem 
when there is an upper bound on the radii. Section 4 provides a number of 
approximation results in case there is no upper bound on the radii. In Section 5 
we present a PTAS for the general case, complemented by experiments in Section 
6. A concluding discussion with open problems is provided in Section 7. 

2 CRA for a Given Connectivity Tree 

For a given connectivity tree, our problem is polynomially solvable, based on the 
following observation. 

Lemma 1. Given a connectivity tree T with at least three nodes. There exists 
an optimal range assignment for T with = for all leaves Pi of T . 

Proof. Assume an optimal range assignment for T has a leaf pi € P with radius 
r.i > 0. The circle C; around pi with radius intersects circle Cj around p^s 
parent pj with radius r,-. Extending Cj to rj := dist(j>i,pj) while setting r*j := 
does not increase the solution value R — X) P ep r i- ^ 

Direct consequences of Lemma [l] are the following. 

Corollary 1. There is an optimal range assignment satisfying Lemma^ and 
rj > for all pj € P of height 1 in T (i.e., each pj is a parent of leaves only). 

Corollary 2. Consider an optimal range assignment for T satisfying Lemma^ 
Further let pj £ P be of height 1 inT. Thenrj > max Pi , s c hud of pj{dist(pi,pj)}. 

These observations allow a solution by dynamic programming. The idea is to 
compute the values for subtrees, starting from the leaves. Details are omitted. 

Theorem 1. For a given connectivity tree, CRA is solvable in 0(n). 

3 Range Assignment for Bounded Radii 

Without a connectivity tree, and assuming an upper bound of p on the radii, 
the problem becomes NP-hard. In this extended abstract, we sketch a proof of 
NP-hardness for the graph version of the problem; for the geometric version, a 
suitable embedding (based on Planar 3SAT) can be used. 

Theorem 2. With radii bounded by some constant p, the problem CRA is NP- 
hard in weighted graphs. 

See Figure [l] for the basic construction. The proof uses a reduction from 
3Sat. Variables are represented by closed "loops" at distance p that have two 
feasible connected solutions: auxiliary points ensure that either the odd or the 
even points in a loop get radius p. (In the figure, those are shown as bold black 
or white dots. The additional small dots form equilateral triangles with a pair 
of black and white dots, ensuring that one point of each thick pair needs to be 
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Fig. 1. Two variable gadgets connected to the same clause gadget. "True" and "False" 
vertices marked in bold white or black; auxiliary vertices are indicated by small dots; 
the clause vertex is indicated by a triangle. Connectivity edges are not shown. 



chosen, so a minimum-cardinality choice consists of all black or all white within 
a variable.) Additional "connectivity" edges ensure that all variable gadgets are 
connected. Each clause is represented by a star-shaped set of four points that is 
covered by one circle of radius p from the center point. This circle is connected 
to the rest of the circles, if and only if one of the variable loop circles intersects 
it, which is the case if and only if there is a satisfying variable. 

4 Solutions with a Bounded Number of Circles 

A natural class of solutions arises when only a limited number of k circles may 
have positive radius. In this section we show that these k-circle solutions already 
yield good approximations; we start by giving a class of lower bounds. 

Theorem 3. A best k-circle solution may be off by a factor of (1 + 2fc+ 1 i_ 1 ). 



Radii: 2°2 1 2 2 2 k ~ l , 2 




Fig. 2. A class of CRA instances that need k + 1 circles in an optimal solution. 



Proof. Consider the example in Fig. [2] The provided solution r is optimal, as 
R :— Ti — di st (Po-P") _ Further, for any integer k > 2 we have d\ — 2 J + 

2 fe ~ 1 < 2 • 2 k + 2 fc ~ 1 = d 2 . So the radius r^ + i cannot be changed in an optimal 
solution. Inductively, we conclude that exactly k + 1 circles are needed. Because 
we only consider integer distances, a best fc-circle solution has cost Rk > R+ 1, 
i.e., f> 1+2^31. □ 
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In the following we give some good approximation guarantees for CRA using 
one or two circles. 

Lemma 2. Let V a longest (simple) path in an optimal connectivity graph, and 
lete m be an edge inV containing the midpoint ofV. Then^ri > max{||7-'|, |e m |} 

This follows directly from the definition of the connectivity graph which for any 
edge e = p u Pv in V requires r u + r v > \e\. 

Theorem 4. A best l-circle solution for CRA is a ^-approximation, even in 
the graph version of the problem. 

Proof. Consider a longest path'P = (po, ■ ■ ■ ,Pk) of length \V\ = d-p(po, ■ ■ ■ ,Pk) '■= 
Yli^o \PiPi+i I i n the connectivity graph of an optimal solution. Let R* :— Yl r * 
be the cost of the optimal solution, and e m = piPi+i as in Lemma[2] Let di := 
d-pipu ■ ■ ■ ,Pk) and d i+1 := d v (p , . . . ,JJ»+i). Then min{Ji,J i+ i} < ' + 2 ' +1 = 
d y ( P o,-,p i )+2|^|+<MPi + i,-,p fc ) = M + < R * + B*. = 3 R *. So one circle with 
radius around the point in P that is nearest to the middle of path V covers 
P, as otherwise there would be a longer path. □ 



optimal solution 
best l-circle solution 

Fig. 3. A lower bound of § for l-circle solutions. 

Fig. [3] shows that this bound is tight. Using two circles yields an even better 
approximation factor. 

Theorem 5. A best 2-circle solution for CRA is a ^-approximation, even in 
the graph version of the problem. 

Proof. Let V = (p , . . . ,pk) be a longest path of length \V\ = d-p(p , ■ ■ ■ ,Pk) '■= 
X^c) 1 |PiP*+i| m the connectivity graph of an optimal solution with radii r*. 
Then R* :=J2 r i ^ ll^l- We distinguish two cases; see Fig. [i] 

Case 1. There is a point x on V at a distance of at least \\V\ from both 
endpoints. Then there is a l-circle solution that is a |-approximation. 

Case 2. There is no such point x. Then two circles are needed. One of 
them is placed at a point in the first third of V , and the other circle is placed 
at a point in the last third of V . Let e m = piPi+\ be defined as in Lemma [2j 
Further, let di := d v (po, ■ ■ ■ ,Pi), and let d l+i := d-p(p i+1 , . . . ,p k ). Then \e m \ = 
\V\ - di - d l+ i and d { ,d i+1 < \\V\. 
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Fig. 4. The two |-approximate 2-circle solutions constructed in the proof of Theorem[5j 
(Top) case 2a; (bottom) case 2b. 



Case 2a. If |e m | < \\V\ then d t + d l+1 = \V\ - \e m \ > \\V\ > \e m \. Set 
Ti :— di and r^+x := di+i, then the path is covered. Since di,di + \ < ^{Vl we 
have Ti + r,+x = <i; + di + i < ^\V\ < f-R* and the claim holds. 

Case 2b. Otherwise, if |e m | > ^\V\ then d{ + d i+ i < \\V\ < \e m \. Assume 
di > di+%. Choose := di and rj+i := |e m | — d^. As di + \ < \e m \ — di the path 
V is covered and r.; + r,-+x = di + (\e m \ — di) = |e m |, which is the lower bound 
and thus the range assignment is optimal. □ 

If all points of P lie on a straight line, the approximation ratio for two circles 
can be improved. 

Lemma 3. Let P be a subset of a straight line. Then there is a non-overlapping 
optimal solution, i.e., one in which all circles have disjoint interior. 

Proof. An arbitrary optimal solution is modified as follows. For every two over- 
lapping circles Ci and Cj+x with centers pi and Pi+\, we decrease r^+x, such that 
r i + r i+i = dist(pi,Pi + i), and increase the radius of Ci + 2 by the same amount. 
This can be iterated, until there is at most one overlap at the outermost circle Cj 
(with Cj-i). Then there must be a point Pj+\ on the boundary of Cj\ otherwise 
we could shrink Cj contradicting optimality. Decreasing Cj's radius Vj by the 
overlap I and adding a new circle with radius I around Pj+\ creates an optimal 
solution without overlap. □ 

Theorem 6. Let P a subset of a straight line g. Then a best 2-circle solution 
for CRA is a ^-approximation. 

Proof. According to Lemma [3] we are, w.l.o.g., given an optimal solution with 
non-overlapping circles. Let po and p n be the outermost intersection points of the 



Connecting a Set of Circles with Minimum Sum of Radii 

R* , R* 



+ 



H 
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optimal solution circles and g. W.l.o.g., we may further assume Po,p n € P and 



R* 



dist(p 



(otherwise, we can add the outermost intersection point 



of the outermost circle and g to P, which may only improve the approximation 
ratio) . Let Pi denote the rightmost point in P left to the middle of Pop n and 
let Pi+i its neighbor on the other half. Further, let di : = dist(pojPi)) di+i '■= 
dist(pi+i,p n ) (See Fig. [5]). Assume, di > dj+i. We now give |-approximate 
solutions using one or two circles that cover pop n - 

Case 1. If fiT < di then \R* > 2R* - d % = dist(p;,p„). Thus, the solution 
consisting of exactly one circle with radius 2R* — di centered at pi is sufficient. 

Case 2. If %R* > di > d^+i we need two circles to cover pop n with \R*. 



Po\ 



Pi 



1 r>* 1 D* 



Pi Pr 



di ' \R* — di R* — di 



Fig. 6. A |-approximate 2-circle solution with di < f -R*. The cross marks the position 
of the optimal counterpart p* to pi and the grey area sketches Ai. 



Case 2a.Thc point pi could be a center point of an optimal two-circle solution 
if there was a point p* with dist((7j,p*) = dist(p*,p„) = R* — d t . So in case 
there is a p\ e P that lies in a -neighborhood of such an optimal p* we get 
dist((7 i ,p' j ) ) dist(pi,p n ) < R* - d t + \R* (see Fig.g. Thus, r{ Pi ) := d t ,r{p' % ) := 
R* — di + jR* provides a | -approximate solution. 



8 Chambers, Fekete, Hoffmann, Marinakis, Mitchell, Srinivasan, Stege, et al. 



Case 2b. Analogously to Case 2a, there is a point p' i+1 G P within a \R*- 
range of an optimal counterpart to Pi+i- Then we can take r(pi+i) :— d^+i, 
r(p' i+1 ) := R* — d i+ i + ^R* as a |-approximate solution. 

Case 2c. Assume that there is neither such a p[ nor such a p' i+1 - Because 
di, d i+ i are in (\R* , jR*), we have \R* < R* - dj < for j = i,i + 1, which 
implies that there are two disjoint areas Ai, Ai+i, each with diameter equal to 
t^R* and excluding all points of P. Because pi, the rightmost point on the left 
half of poPn, has a greater distance to Ai than to po, any circle around a point 
on the left could only cover parts of both Ai and A i+ i if it has a greater radius 
than its distance to pq. This contradicts the assumption that po is a leftmost 
point of a circle in an optimal solution. The same applies to the right-hand side. 
Thus, Ai U Ai + i must contain at least one point of P, and therefore one of the 
previous cases leads to a |-approximation. □ 




Fig. [7] shows that the bound is tight. We believe that this is also the worst 
case when points are not on a line. Indeed, the solutions constructed in the proof 
of Theorem [6] cover a longest path V in an optimal solution for a general P. If 
this longest path consists of at most three edges, Pi{=: p'i + \) and pi + \{—: p'i) can 
be chosen as circle centers, covering all of P. However, if V consists of at least 
four edges, a solution for the diameter may produce two internal non-adjacent 
center points that do not necessarily cover all of P. 

5 Polynomial-Time Approximation Schemes 

We now consider the problem in which each of the n points of P = {pi, . . . 
has an associated upper bound, fi, on the radius Ti that can be assigned to pi. 

5.1 Unbounded Radii 

We begin with the case in which fj = oo, for each i. Consider an optimal solution, 
with radius r* associated with input point p t . We first prove a structure theorem 
that allows us to apply the m-guillotine method to obtain a PTAS. The following 
simple lemma shows that we can round up the radii of an optimal solution, at a 
small cost to the objective function: 



Connecting a Set of Circles with Minimum Sum of Radii 



9 



Lemma 4. Let R* — r* be the sum of radii in an optimal solution, T>* . 
Then, for any fixed e > 0, there exists a set, T> m , of n circles of radii j-j centered 
on points Pi, such that (a), r, € K = {D /mn,2D /mn, . . . ,D}, where D is the 
diameter of the input point set P and m — |~2/e~|; and (b). J^. r\ < (l + e)R*. 

Proof. Each of the n radii r* can be increased by at most D /mn < eD/2n at a 
total cost of at most eD/2. Since increasing the radii of the circles keeps the set 
of circles connected and since R* > D/2, we obtain the result. □ 

Disks centered at the points P of radii in the set 1Z = {D/mn, 2D/mn, . . . , D} 
will be referred to as lZ e ,p -circles, or IZ-circles, for short, with the understand- 
ing that e and P will be fixed throughout our discussion. Consider the arrange- 
ment of all TvL-circles. We let l x (resp., I y ) denote the ^-coordinates of the 
left/right (resp., ^/-coordinates of the top/bottom) extreme points of these cir- 
cles. (Specifically, I x = {x Pi ±j(D/mn) : 1 < i < n, < j < mnfi/D} and 
X y = {y p . ±j(D/mn) : 1 < i < n,0 < j < mnfi/D}.) 

We say that a set T> of n 7?.-circles is m-guillotine if the bounding box, 
BB(D), of T> can be recursively partitioned into a rectangular subdivision by 
axis-parallel "m-perfect cuts" that are defined by coordinates I x and X y , with the 
finest subdivision consisting of a partition into rectangular faces each of which 
has no circle of T> strictly interior to it. An axis-parallel cut line i is m-perfect 
with respect to V and a rectangle p if £ intersects at most 2m circles of T> that 
have a nonempty intersection with p. 

Key to our method is a structure theorem, which shows that we can transform 
an arbitrary set T> of circles centered on points P, having a connected union and 
a sum of radii R, into an m-guillotine set of 7?.-circles, T> m , having sum of radii 
at most (1 + e)R* . More specifically, we show (proof deferred to the full paper): 

Theorem 7. Let T> be a set of circles of radii centered at points pi G P , such 
that the union of the circles is connected. Then, for any fixed e > 0, there exists 
an m-guillotine set T> m of n IZ-circles such that the union of the circles T> m is 
connected and the sum of the radii of circles ofT> m is at most (1 + (C/m)) r^. 
Here, m = fl/e] and C is a constant. 

Proof. First, by Lemma |4j we can afford to round up the radii of the circles T> 
to make them 7?.-circles. 

Now, let p denote the axis-aligned bounding box of the resulting circles. If 
there are no circles strictly interior to p (i.e., if all circles intersecting p intersect 
(touch) the boundary, dp, of p), then we are done: the set of circles is trivially m- 
guillotine already. Thus, we assume that there is at least one circle of T> strictly 
interior to p. 

If there exists an m-perfect cut of p, by a horizontal or vertical line inter- 
secting at most 2m circles within p and intersecting at least one circle that lies 
interior to p, then we partition p with it and recurse on the two boxes on each 
side of the cut. We observe that if such an m-perfect cut exists, then one exists 
with the property that the defining coordinate of the cut is from the set I x (for 
vertical cuts) or I y (for horizontal cuts), since we can translate the cut between 
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two consecutive such coordinates without changing the set of circles it intersects. 
(More precisely, we can translate a cut to be infinitesimally close to one of the 
discrete coordinates without changing its combinatorial type. To address this 
technicality, we can either augment the sets T x and T y with the midpoints of 
the intervals between consecutive coordinates, or we can replace each coordinate 
z with three coordinates - z and z + (infinitesimally greater than z) and z~ 
(infinitesimally smaller than z).) 

Thus, we now assume that no m-perfect cut exists for partitioning p, and 
that p contains in its interior at least one circle. For a vertical line, £ x , through 
coordinate x, let f(x) denote the length of the m-span of £ x with respect to V 
and p: f(x) = if, within p, £ x intersects at most 2m circles of V\ otherwise, if 
l x intersects K > 2m circles of T> within p, then f(x) is the distance (along £ x ) 
from the first point, a m , where £ x enters into the mth circle, going from the top 
boundary of p downwards along £ x , to the first point, b m , where £ x enters into 
the mth circle, going from the bottom boundary of p upwards along £ x . Because 
K > 2m, we know that a m is above b m and that there must be at least one circle 
of T> whose intersection with £ x is a proper subset of the bridge segment a m b m . 
We similarly define g(y) to be the length of the (horizontal) bridge segment along 
a horizontal cut £ y through coordinate y. 

We think of f(x) and g(y) as the cost of augmentation for the network JV p 
consisting of the union of circles, truncated within p, that are the boundaries 
of the circles V; by adding segments (bridges) of length f(x) (resp., g(yj), a 
vertical cut £ x (resp., horizontal cut £ y ) can be made m-perfect. We claim that 
we can charge off the lengths of the bridges that would suffice to augment the 
network J\f p to make it m-guillotine, in the usual sense of an m-guillotine network 
(subdivision), as in |11) . Specifically, we argue that we can select cuts for which 
the bridge lengths (m-spans) can be charged off to the total length of the network 
(sum of the circle circumferences, which is 0(X^ r i))i showing that the sum of 
the bridge lengths is at most (C/m) ^ r^. 

We partition each circle (bounding the circles V) into four 90-degree arcs: two 
"vertical arcs" (with angular ranges (-45,45) and (135,225)) and two "horizontal 
arcs" (with angular ranges (45,135) and (225,315)). We define the "chargeable 
length" of a vertical cut £ x to be the "m-dark" length of £ x H p. Specifically, a 
subsegment ab of £ x n p is said to be m-dark with respect to Af if for any p £ a&, 
the rightwards and leftwards rays from p each cross at least m vertical arcs of 
Af p before exiting p. If we cut p along £ x , then the m-dark portion of the cut 
can be charged off to the left /right sides of the vertical arcs of Af p lying to the 
right/left of £ x , distributing the charge to be (l/m)th to each of the m arcs first 
hit. 

Our charging scheme is based on the observation, following the method of 
pTj . that there must exist a "favorable" vertical cut £ x or horizontal cut £ y 
for p such that the chargeable length of the cut is at least as long as the cost 
(f(x) or g(y)) of the cut. The existence of a favorable cut follows from the ob- 
servation that J x( - p f(x)dx = f h(y)dy, where h(y) is the chargeable length 
associated with the horizontal cut £ y ; thus, assuming, without loss of generality, 
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that J xep f(x)dx > J yep g(y)dy, we see that there must exist a value y* where 
g(y*) < h{y*), which defines a favorable horizontal cut l y * for which the charge- 
able length exceeds the length of the m-span. (In case J xep f(x)dx < J ye g(y)dy, 
there exists a favorable vertical cut.) Further, we claim that there must exist a 
favorable cut corresponding to the coordinate sets l x ,l y : the chargeable length 
of cut ly- does not change as we perturb y* to the nearest coordinate in the set 
l y , while, by convexity of the circular arcs, the length of the m-span associated 
with a cut is locally minimized at endpoints of the intervals defined by the points 

Of ly. 

Once a favorable cut is found with respect to rectangle p, the cut partitions 
the problem into two subrectangles, and the argument is recursively applied to 
each. Since each circular arc of length A is charged for length at most A/2m 
on each of its two sides, we get that the overall length of all m-spans that are 
associated with favorable cuts constructed recursively in converting the network 
Nbb(d) to an m-guillotine network is at most (C/m) J2i r ii f° r a constant C. 

Finally, we claim that the circles T> can be converted to an m-guillotinc 
set T> m , having sum of radii at most (C /m) Y], Tj : Associated with each m- 
span a m b m that is added to the network of circular arcs in order to make the 
network m-guillotine, we enlarge the radius of the circle defining one of the m- 
span endpoints (say, a m ), by at most |a m & m |, so that this circle now covers the 
entire m-span segment. Since we have enlarged one circle in a set of circles whose 
union is connected, the union remains connected. Thus, with a total increase in 
radii of at most (C/m) rj, we end up with an m-guillotine set T> m of circles, 
proving our structure theorem. □ 

We now give an algorithm to compute a minimum-cost (sum of radii) m- 
guillotine set of 7£-circles whose union is connected. The algorithm is based on 
dynamic programming. A subproblem is specified by a rectangle, p, with x- and 
y-coordinates among the sets l x and l y , respectively of discrete coordinates. 
The subproblem includes specification of boundary information, for each of the 
four sides of p. Specifically, the boundary information includes: (i) O(m) "portal 
circles", which are 7?.-circles intersecting the boundary, dp, of p, with at most 
2m circles specified per side of p; and, (ii) a connection pattern, specifying which 
subsets of the portal circles are required to be connected within p. There are a 
polynomial number of subproblems, for any fixed m. For a given subproblem, 
the dynamic program optimizes over all (polynomial number of) possible cuts 
(horizontal at ^-coordinates or vertical at X„-coordinates) , and choices of up to 
2m "/^-circles intersecting the cut bridge, along with all possible compatible con- 
nection patterns for each side of the cut. The result is an optimal m-guillotine 
set of 7?.-circles such that their union is connected and the sum of the radii 
is minimum possible for m-guillotine sets of 7^-circles. Since we know, from the 
structure theorem, that an optimal set of circles centered at points P can be con- 
verted into an m-guillotine set of 7?.-circles centered at points of P, whose union 
is connected, and we have computed an optimal such structure, we know that 
the circles obtained by our dynamic programming algorithm yield an approxi- 
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mation to an optimal set of circles. In summary, we have shown the following 
result: 

Theorem 8. There is a PTAS for the min-sum radius connected circle problem 
with unbounded circle radii. 

5.2 Bounded Radii 

We now address the case of bounded radii, in which circle i has a maximum 
allowable radius, fj < oo. The PTAS given above relied on circle radii being 
arbitrarily large, so that we could increase the radius of a single circle to cover the 
entire m-span segment. A different argument is needed for the case of bounded 
radii. 

We obtain a PTAS for the bounded radius case, if we make an additional 
assumption: that for any segment pq there exists a connected set of circles, 
centered at points of pi € P and having radii 7"j < f,-, such that p and q each lie 
within the union of the circles and the sum of the radii of the circles is 0(|pg|). 

Here, we only give a sketch of the method, indicating how it differs from the 
unbounded radius case. The PTAS method proceeds as above in the unbounded 
radius case, except that we now modify the proof of the structure theorem by 
replacing each m-span bridge a m b m by a shortest connected path of TvL-circles. 
We know, from our additional assumption, that the sum of the radii along such 
a shortest path is 0(\a m b m \), allowing the charging scheme to proceed as before. 
The dynamic programming algorithm changes some as well, since now the sub- 
problem specification must include the "bridging circle-path" , which is specified 
only by its first and last circle (those associated with the bridge endpoints a m 
and b m ); the path itself, which may have complexity fi(n), is implicitly specified, 
since it is the shortest path (which we can assume to be unique, since we can 
specify a lexicographic rule to break ties). 

In summary, we have 

Theorem 9. There is a PTAS for the min-sum radius connected circle problem 
with bounded circle radii, assuming that for any segment pq, with p and q within 
feasible circles, there exists a (connected) path of feasible circles whose radii are 
0(\ P q\). 

6 Experimental Results 

It is curious that even in the worst case, a one-circle solution is close to be- 
ing optimal. This is supported by experimental evidence. In order to generate 
random problem instances, we considered different numbers of points uniformly 
distributed in a 2D circular region. For each trial considering a single distri- 
bution of points, we enumerated all possible spanning trees using the method 
described in [2], and recorded the optimal value with the algorithm mentioned 
in Section [2] This we compared with the best one-circle solution; as shown in 
Fig. [HJ the latter seems to be an excellent heuristic choice. These results were 
obtained in several hours using an i7 PC. 
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Fig. 8. Ratios of the average over all enumerated trees and of the best 1-circle tree to 
the optimal ^2 ri. Results were averaged over 100 trials for each number. 

7 Conclusion 

A number of open problems remain. One of the most puzzling is the issue of 
complexity in the absence of upper bounds on the radii. The strong performance 
of the one-circle solution (and even better of solutions with higher, but limited 
numbers of circles), and the difficulty of constructing solutions for which the 
one-circle solution is not optimal strongly hint at the possibility of the problem 
being polynomially solvable. Another indication is that our positive results for 
one or two circles only needed triangle inequality, i.e., they did not explicitly 
make use of geometry. 

One possible way may be to use methods from linear programming: modeling 
the objective function and the variables by linear methods is straightforward; 
describing the connectivity of a spanning tree by linear cut constraints is also 
well known. However, even though separating over the exponentially many cut 
constraints is polynomially solvable (and hence optimizing over the resulting 
polytopc), the overall polytope is not necessarily integral. On the other hand, 
we have been unable to prove NP-hardness without upper bounds on the radii, 
even in the more controlled context of graph-induced distances. Note that some 
results were obtained by means of linear programming: the tight lower bound 
for 2-circle solutions (shown in Fig. 7) was found by solving appropriate LPs. 

Other open problems are concerned with the worst-case performance of heuris- 
tics using a bounded number of circles. We showed that two circles suffice for a 
|-approximation in general, and a |-approximation on a line; we conjecture that 
the general performance guarantee can be improved to | , matching the existing 
lower bound. Obviously, the same can be studied for k circles, for any fixed k; at 
this point, the best lower bounds we have are | for k — 3 and 1 + for general 
k. We also conjecture that the worst-case ratio f(k) of a best /c-circle solution 
approximates the optimal value arbitrarily well for large k, i.e., lim^oo f(k) = 1. 
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